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ABSTRACT 

The  rf-excited  coil  in  the  presence  of  a  plasma  is  examined  in  terms  of  a 
transformer  model  and  also  in  terms  of  the  distortion  of  the  electromagnetic  fields  which 
is  produced  by  the  induced  eddy-currents  in  the  plasma.  The  treatment  of  rf  coils 
presented  here  can  be  applied  to  either  the  plasma  diagnostics  regime  or  to  the  generation 
of  the  rf  electrodeless  discharge.  Although  the  calculations  of  the  electric  and  magnetic 
fields  in  this  report  are  for  a  cylindrically  symmetric  coil  immersed  in  a  uniform,  linear 
plasma,  neither  the  model  nor  the  approach  to  the  problem  is  limited  to  these  cases. 
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V  Potential  across  the  coil 
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SECTION  I 
INTRODUCTION 

Interest  in  the  coupling  between  radio-frequency  (if)  excited  coils  and  conductive 
media  stems  from  the  fact  that  rf-exdted  coils  are  used  both  for  generation  of  plasmas 
and  for  plasma  diagnostics.  In  the  diagnostic  regime,  rf  coils  have  been  used  to  measure 
plasma  parameters  such  as  electrical  conductivity,  gas  flow  velocity,  and  electron  number 
density  (Refs.  1,  2,  and  3).  RF  induction  heating  of  metals  has  long  been  a  method  of 
heat  treating  metals,  and  recently  gas  plasmas  have  been  heated  to  high  enthalpy  levels 
using  rf  induction  techniques  (Refs.  4,  5,  and  6).  Although  there  is  a  gulf  which  separates 
these  two  areas  of  application  of  rf  coils  to  plasma  physics,  there  is  a  common 
denominator  which  links  the  two  phenomena,  namely,  both  are  induced  effects  created 
by  a  time  varying  magnetic  field.  In  the  former  case  (i.e.,  plasma  diagnostics),  the  fields 
are  very  weak,  and  the  distortion  of  these  fields  by  the  presence  of  the  plasma  is  related 
to  various  properties  of  the  plasma.  In  the  latter  case,  the  fields  are  much  stronger,  and 
the  coupling  mechanism  is  much  more  involved  since  the  conductivity  of  the  plasma  is 
not  uniform.  However,  a  great  deal  of  useful  insight  can  be  obtained  by  treating  the 
system  as  though  it  were  uniform  and  isotropic. 

Since,  in  the  final  analysis,  it  becomes  necessary  to  explain  the  interaction 
between  the  coil  and  the  plasma  in  terms  of  circuit  phenomena  (e.g.,  current,  potential, 
resistance,  inductance,  and  capacitance),  it  is  expedient  to  employ  from  the  beginning  a 
model  which  characterizes  the  problem  in  terms  of  these  quantities.  It  should  be  pointed 
out  also  that  the  model  to  be  presented  here  does  not  require  that  the  conductivity  of 
the  plasma  be  uniform.  Rather,  this  requirement  is  used  only  to  make  some  simple 
analytical  calculations  for  illustrational  purposes.  A  host  of  experimental  evidence 
supports  the  use  of  a  coupling  model  in  which  the  plasma-coil  relationship  is 
approximated  by  a  simple  transformer  circuit  (Refs.  7,  8,  and  9).  In  this  model,  the 
distribution  of  the  induced  eddy-currents  is  treated  in  terms  of  lumped  circuit  parameters 
which  can  be  reflected  into  the  primary  circuit  using  ordinary  transformer  techniques. 
Hence,  a  qualitative  analysis  of  the  changes  in  the  coil  inductance  and  resistance  in  the 
presence  of  the  plasma  can  be  made.  Exact  values  of  the  coil  parameters  in  the  presence 
of  a  plasma  is  obtained  by  solution  of  Maxwell’s  equation.  In  some  simple  geometrical 
configuration  and  for  uniform  conductivity  or  for  some  special  cases  (e.g.,  parabolic 
conductivity,  Ref.  6),  analytical  expressions  for  these  configurations  are  readily 
obtainable,  whereas  in  the  more  complex  cases,  a  numerical  solution  must  be  obtained. 

Although  the  solution  to  the  coupling  problem  must  be  sought  by  examining  the 
induced  eddy-currents,  regardless  of  any  particular  model,  the  method  in  which  this 
problem  is  treated  can  be  simplified  by  seeking  a  correlating  parameter  which  can  be 
readily  interpreted  in  terns  of  either  circuit  phenomena  or  field  phenomena.  The  “figure 
of  merit”  or  Q  of  a  circuit  (defined  as  the  ratio  of  average  energy  stored  to  the  energy 
loss  per  cycle)  is  a  nondimensional  parameter  which  is  ideal  for  characterizing  the 
eddy-currents,  and  it  is  also  a  parameter  which  is  directly  applicable  to  explaining  circuit 
phenomena.  The  coupling  between  rf-ex cited  coils  and  conductive  media  is  explained  in 
this  report  by  exploiting  the  features  associated  with  the  figure  of  merit  of  the  system 
(Qs)  in  conjunction  with  a  simple  transformer  model.  A  comparison  between  theory  and 
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experiment  for  solenoids  immersed  in  conductive  media  indicates  that  the  coupling  can 
be  accurately  explained  (both  qualitatively  and  quantitatively)  using  the  transformer 
model. 


It*  should  be  noted  that  the  analytical  and  experimental  results  presented  in  this 
report  do  not  depend  on  the  transformer  model.  This  model  is  intended  primarily  for 
clarification  and  understanding.  Similarly,  the  transformer  model  is  not  restricted  to  the 
simple  case  to  be  studied  here. 


SECTION  II 

TRANSFORMER  MODEL 

It  has  been  observed  that,  when  a  nonmagnetic  conductor  is  placed  in  the 
neighborhood  of  an  rf-exdted  inductor,  (1)  eddy-currents  are  induced  in  the  conductor, 
(2)  the  self  inductance  of  the  inductor  decreases,  and  (3)  the  internal  resistance 
associated  with  the  inductor  increases.  These  effects  are  also  produced  when  a  closed 
secondary  circuit  is  inductively  coupled  to  the  inductor.  This  parallelism  suggests  a  simple 
model  which  can  be  used  to  (qualitatively)  explain  the  coupling  mechanism  between  rf 
coils  and  conductive  media  (plasmas).  That  is,  the  eddy-currents  may  be  treated  by 
lumping  the  distributed  inductance  and  resistance  associated  with  each  closed  current 
path  into  a  total  equivalent  inductance  (L|)  and  total  equivalent  resistance  (R*).  Further, 
if  it  is  assumed  that  the  coil  is  inductively  coupled  to  the  plasma  through  a  coupling 
coefficient  k,  then  the  effects  caused  by  the  presence  of  the  plasma  can  be  easily 
examined  by  reflecting  the  lumped  inductance  and  resistance  into  the  coil  using  ordinary 
transformer  analysis.  This  procedure  is  schematically  illustrated  in  Fig.  1 . 


L  * 

R  * 


Fig.  1  Equivalent  Circuit  Model  of  an  RF-Excited  Coil  in  the  Presence 
of  a  Plasma 


From  this  figure,  it  can  be  seen  that  the  process  is  divided  into  three  successive 
steps.  First,  the  coil  with  free-space  inductance  L  and  resistance  R  is  shown  immersed  in 
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a  conductive  medium.  In  the  second  step,  the  inductance  and  resistance  associated  with 
the  induced  eddy-currents  are  treated  as  a  closed  secondary  circuit  coupled  to  the 
primary  circuit  through  a  coupling  coefficient  k.  In  the  third  step,  the  secondary  circuit 
is  reflected  into  the  primary  circuit  which  results  in  a  new  coil  whose  inductance  and 
resistance  (1/  and  R')  can  be  obtained  directly  from  the  loop  equations  for  the  circuit 
shown  in  the  center  of  Fig.  1.  That  is, 

V  =  (R  +  jfijL)  Ip  +  {jiukVLLB)  Ia  (1) 

0  =  (jii>k.\/LLs)  Ip  +  (Rs  +■  jct)Ls)  Ift  (2) 

where  V  is  the  applied  potential  across  points  (1)  and  (2)  in  Fig.  1,  or  rearranging  Eq.  (2) 
gives 


Substituting  for  Ig  in  Eq.  (3)  and  taking  the  ratio  V/Ip  yield  the  input  impedance  Z;n  of 
the  coil: 


Zin  «  - 


1  + 


Ic2<wL)(<uLb)  R. 


R[R,2+  (<uL  b  ) 2] 


+  jwL 


1  - 


k2  (^L,)2 


R,2  +  (&>LB)2 


Equation  (4)  can  be  written  in  the  form 

Zjn  =  R  +  jwL 


where 


L'  =  LI 


k2(<aL.)2 
R  l  +  (coL,)2 


and 


(4) 

(5) 

(6) 


R'  -  R 


1  + 


k2(cuL)  (<aLs)  Rb 
R[Rf  +  (cjLb)2] 


(7) 


If  the  nondimensional  parameters  Q  =  wL/R,  Q'  =  cuL'/R',  and  Qs  =  ojLs/R,  are 
substituted  into  Eqs.  (6)  and  (7),  then  the  following  useful  relationships  are  obtained: 


L' 

L 


l  +  q2(i  -kz) 
l  +  Q* 


<  1 


(8) 


R'  q2  +  k2Q  Qa  +  1 
R  ~  Q?+1 


(9) 
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and 


51 .  n-oju-t*?  , 
Q  q2  +i5Q  Q.  +  l 


CIO) 


The  figure  of  merit  Q,  which  is  associated  with  the  induced  eddy-currents  in  the  plasma 
can  be  expressed  in  terms  of  the  measurable  quantities  L,  L\  R,  and  R'  by  eliminating 
the  coupling  coefficient  k  in  Eqg.  (8)  and  (9).  That  is, 


Q. 


<u(L  -  L ')  /Al\ 

R'-R  \AR/ 


(ID 


The  average  power  input  into  the  plasma-coil  system  can  be  expressed  in  terms  of 
L,  Q,  Q, ,  k,  and  the  current  Ip : 

Pin  -  Pc  oil  +  Ppl  asm  a  =  IpR '  (12) 

where 


A  number  of  conclusions  regarding  the  application  of  rf-ex cited  coils  to  specific 
problem  areas  can  be  made  from  examination  of  Eq.  (14).  For  example,  it  can  be  seen 
that,  for  diagnostic  applications,  it  is  necessary  that  Pcon  be  small  (of  the  order 
microwatts)  so  that  the  plasma  is  undisturbed,  while  for  rf  plasma  generation,  it  is 
desirable  to  have  PpitBm,  »Pcoii  (i-e.,  Q,  =  1  and  k2Q»  1)  and  Ppitlffl»  to  be  of  the 
order  of  kilo-  to  megawatts.  In  addition  to  these'  considerations,  it  is  also  important  to 
consider  matching  between  the  coil  and  electronics  as  well  as  matching  between  the  coil 
and  the  plasma.  A  brief  discussion  of  this  matching  is  presented  in  Appendix  I. 

SECTION  lit 

ANALYTICAL  PROCEDURES 


3.1  FIGURE  OF  MERIT  Q, 

The  emphasis  to  this  point  has  been  to  describe  the  coupling  between  the  coil  and 
plasma  in  terms  of  circuit  parameters.  This  is  useful  because,  in  the  final  analysis, 
measurements  of  the  system  are  greatly  simplified  if  everything  can  be  put  in  terms  of 
linear  circuit  elements.  However,  while  this  approach  is  desirable  in  establishing  methods 
of  measuring  electrical  conductivity  or  in  measuring  such  things  as  the  input  power  to 
rf-geneiated  plasmas,  it  does  not  lead  to  a  direct  relationship  between  the  measured 
quantities  (e.g.,  Ip ,  Q\  L',  or  R')  and  the  electrical  properties  of  the  plasma  (e.g.,  a,  e, 
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etc.).  This  portion  of  the  report  will  be  concerned  with  calculating  the  coil  parameters  in 
terms  of  the  geometry  and  plasma  properties.  First,  a  cylindrical  coil  will  be  treated  in 
which  the  plasma  is  external  to  the  solenoid.  Then,  the  same  analysis  will  be  applied 
where  the  plasma  is  in  the  interior  of  the  coil. 

The  relationship  between  circuit  phenomena  and  field  phenomena  can  be 
established  by  expressing  Qs  in  terms  of  the  electromagnetic  fields  generated  by-  the 
eddy-currents.  The  defining  equation  is 


Time  averaged  energy  stored 

=  4  rr 

KL.i* 

Wd 

Energy  loss  per  cycle 

2rr  o 

- H.  I, 

_  CO  _ 

3.2  AZIMUTHAL  AND  AXIAL  ELECTRIC  FIELDS  ASSOCIATED 
WITH  A  SOLENOID 

The  energy  stored  and  energy  loss  per  cycle  can  be  calculated  when  the  electric 
and  magnetic  fields  are  known  throughout  the  volume  of  interest.  In  fact,  complete 
knowledge  of  the  electric  and  magnetic  fields  is  sufficient  to  determine  all  electric  and 
magnetic  phenomena  which  may  occur  within  this  volume  of  interest.  The  following 
sections  will  be  concerned  with  calculating  these  fields  from  Maxwell’s  equations  for 
certain  prescribed  conditions.  Obviously,  some  major  simplifications  and  limitations  are 
necessary  to  obtain  closed-form  solutions  to  Maxwell’s  equations  for  these  particular 
cases.  For  the  main  part,  these  amplifications  and  limitations  are  noted  during  the 
calculations.  However,  there  is  one  subtle  simplification  which  requires  additional 
consideration  although  the  introduction  of  this  simplification  does  not  alter  the 
calculations  in  this  report.  In  the  calculations  throughout  this  report,  it  is  assumed  that 
the  electric  field  is  derivable  from  only  a  time-varying  magnetic  field  (induction).  (That 
is,  it  is  assumed  that  V*E  =  0.)  Although  this  is  generally  true  for  the  field  space,  it  does 
not  apply  to  the  surface  of  the  coil  since  there  are  time-varying  charges  along  the  coil 
which  give  rise  to  an  electric  field  (potential).  This  field  has  principally  an  axial 
component,  which  is  often  eliminated  by  constructing  a  Faraday  cage1  about  the  coil. 
Since  the  source  impedance  of  this  electric  field  is  quite  large,  these  fields  cannot  deliver 
appreciable  power  into  a  highly  conducting  medium  because  they  are  effectively 
“shorted-out”  by  the  presence  of  a  good  conductor.  This  is  precisely  why  a  Faraday  cage 
can  be  successfully  used  to  eliminate  the  axial  electric  field,  and  also  explains  how  they 
are  eliminated  in  a  high-power  (sometimes  referred  to  as  “inductive”)  rf  discharge 
without  the  addition  of  a  Faraday  cage.  On  the  other  hand,  the  azimuthal  field  is 
associated  with  a  system  which  has  a  low  source  impedance,  and  any  attempt  to  “short” 
these  fields  leads  ultimately  to  rf  induction  heating  since  it  is  possible  to  deliver  large 
amounts  of  power  from  the  load  coil  to  the  conductor  in  this  mode  of  operation.  In  rf 
plasma  discharges,  the  discharge  which  is  the  result  of  the  axial  electric  field  is  often 
referred  to  as  “low-power,”  “capacitive,”  or  “E-field”  discharge,  whereas  the  discharge 
which  is  the  result  of  the  azimuthal  electric  field  is  referred  to  as  “high-power,” 


Faraday  cage  can  be  built  by  placing  conducting  wires  around  or  inside  a  coil  such  that  they  are  in 
alignment  with  the  axial  electric  field  but  arranged  such  that  they  do  not  interfer  with  the  azimuthal  induction  electric 
field. 
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“inductive,”  or  “H-field”  discharge.  The  remainder  of  this  report  will  be  concerned  only 
with  the  “induced  electric”  field  since  the  axial  electric  field  can  either  be  eliminated  or 
neglected. 

3.3  CALCULATION  OF  COIL  PARAMETERS  WHEN  THE  PLASMA  IS 
EXTERNAL  TO  THE  SOLENOID 

To  express  Q,  in  terms  of  the  stored  energy  and  energy  loss,  the  electromagnetic 
fields  associated  with  the  induced  eddy-currents  must  be  calculated  from  Maxwell’s 

equations.  While  it  is  impossible  to  obtain 
a  dosed-form  solution  for  a  general 
current  distribution,  there  are  certain 
geometric  configurations  in  which  the 
symmetry  allows  some  major 
simplifications.  In  particular,  when  an  if 
coil,  whose  length  is  large  compared  with 
its  diameter,  is  immersed  in  a  uniform 
plasma  as  shown  in  Fig.  2,  then  the 
following  statements  are  useful  in 
simplifying  the  problem:  (1)  the  induced 
eddy-currents  have  only  an  azimuthal 
component,  (2)  the  magnetic  field 
produced  by  the  eddy-currents  in  the  free 
space  region  (0<p)  is  uniform,  and  (3) 
the  electric  and  magnetic  fields  in  the 
plasma  vary  predominantly  in  the  radial 
coordinate;  thus  all  derivatives  in  the 
other  coordinates  may  be  neglected. 

Since  the  magnetic  field  produced  by  the  eddy-currents  in  the  region  (0  <  p  <a) 
is  assumed  to  be  constant,  then  the  time-averaged  energy  stored  can  be  calculated  quite 
amply.  That  is, 


z 


Rfl.  2  Two-Dimensional  View  of  e  Solenoid 
Immersed  in  a  Plasma 


where  Bt  is  the  magnetic  field  at  the  boundary  which  is  produced  by  the  eddy  currents, 
and  T  -  2n7w  is  the  period  of  oscillation,  (MKS  units  will  be  used  exclusively  in  this 
report.) 


To  illustrate  that  this  expression  for  W,  is  consistent  with  the  circuit 
interpretation  of  W„  the  integral  form  of  Ampere’s  law  can  be  used  to  obtain  the 
relationship  between  the  secondary  current  I,  (see  Fig.  1)  and  the  magnetic  field. 
Ampere’s  law  is 

^  B.d?  -  J  f  J-nda  -  nDis(t)  (17) 

C  ■ 
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If  the  path  of  integration  C  is  around  the  loop  as  indicated  by  the  points  1,  2,  3,  4,  1, 
in  Fig.  2,  then 

$  B-d£  =  £Ba  =  *>is(0  08) 

1,2, 3, 4,1 


Equation  (18)  can  be  written  in  terms  of  the  rms  current  Is;  that  is, 

VTc. 

Thus,  the  expression  for  the  time  average  of  the  stored  energy  Ws  becomes 


-  fra2*  1  (n 

TO  -  _  =  {  _ 

£ 


wa  = 


a2p<>\  .  1 


T*  _  __  T  T2 
1s  —  p  r-'a  is 


where 


„2  „ 

T  n  a  fl0 

Ls  "  £ 


(19) 


(20) 

(21) 


From  this  expression,  it  is  evident  that  the  equivalent  inductance  associated  with 
the  induced  eddy-currents  in  the  plasma  is  to  a  first  approximation  determined  by  the 
geometry  of  the  system  and  is  independent  of  the  conductivity  of  the  plasma.  Thus,  from 
this  analysis,  it  is  seen  that  the  equivalent  inductance  of  the  plasma  L„  behaves  quite  like 
that  of  ordinary  inductors. 


To  obtain  the  energy  loss  per  cycle,  the  fields  in  the  region  ( p  >  a)  must  be 
obtained  from  Maxwell’s  equations.  The  equations  can  be  simplified  by  making  use  of  the 
fact  that  the  spatial  variation  of  the  fields  normal  to  the  conductive  surface  is  much  more 
rapid  than  variations  parallel  to  the  surface;  hence,  all  derivatives  with  respect  to 
coordinates  parallel  to  the  surface  can  be  neglected  as  compared  with  the  normal 
derivatives.  Thus,  the  gradient  operator  V  may  be  replaced  by  (Ref.  1 0): 


Therefore,  Maxwell’s  equations  for  the  Fields  in  the  plasma  can  be  written  as 


(22) 


and 


(23) 


‘P 


(24) 


In  addition  to  the  assumptions  listed  above,  it  has  been  implicitly  assumed  that  (I)  the 
time  variance  of  the  fields  is  sinusoidal  (i.e.,  e^1),  (2)  the  eddy-current  is  related  to 
the  electric  field  Iic  according  to  Ohm’s  Law  (i.e.,  X  =  cdic),  and  (3)  the  displacement 
current  in  the  conductive  media  can  be  neglected  (i.e.,  a»  cce). 
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The  magnetic  and  electric  fields  in  the  plasma  which  satisfy  Eqs.  (23)  and  (24) 
along  with  the  boundary  conditions  Bc  B  Bt  at  p  =  a  are 


and 


where 

and 


Bc  =  i,  B,  exp  [-a  (/>-■)} 

c-  * .  B,  a  a  ff2  r  ,  u 
Ec  -  - - -  exp  L~a  (p  -  a)J 

«  =  (i  -p/a 

8  =  (2/<u/Ao o)*  -  “■kin-depth” 


(25) 

(26) 


The  energy  dissipated  per  cycle  is 

i/2  itr  « 

Wd  -  (f)}  J  J  J  *  Ec.Ec*p<Wdz 

V  '  i/2  a  • 

Substituting  for  Ec  in  Eq.  (27)  yields 


(27) 


(28) 


If  Eq.  (19)  is  substituted  into  Eq.  (28)  and  then  solved  for  R„  the  following  expression 
for  R,  is  obtained: 


».  - 


2ff  (a  +  S/2) 
ci8 


(29) 


From  this  equation,  it  is  noted  that  the  equivalent  resistance  associated  with  the  induced 
currents  in  the  plasma  is  amply  the  same  resistance  that  is  observed  for  a  uniform 
conductor  whose  cross-sectional  area  is  S6  and  conductance  length  is  2?r(a+  6/2). 

The  figure  of  merit  Q,  for  the  cylindrically  symmetric  system  shown  in  Fig.  2  is 
obtained  by  substituting  Eqs.  (20)  and  (28)  into  Eq.  (15).  That  is, 


where 


Qb 


2 


\v8i_ 

1  +  2b/S 


l  +VT* 


(30) 


x  =  V?  “  aV'fcjpoC 

Equation  (30)  is  plotted  in  Fig,  3  versus  the  nondimensional  parameter  x.  This  figure  also 
contains  some  experimental  data  which  were  obtained  using  a  number  of  different 
cylindrical  coils  immersed  in  various  conductors  whose  conductivities  extended  from  10 
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mhos/m  (electrolytic  solutions)  to  106  mhos/m  (mercury).  The  coil  inductance  and  coil 
Q  were  measured  with  a  Q-metcr  when  the  coil  was  immersed  in  these  conductors  while 
the  frequency  was  varied  between  2  and  20  mHz.  This  measurement  process  is  outlined 
in  detail  in  Ref.  7.  A  comparison  between  theory  and  experiment  indicates  that  the 
transformer  model  adequately  describes  (both  qualitatively  and  quantitatively)  rf  coils 
immersed  in  conductive  media  over  a  wide  range  of  conductivities  and  frequencies. 

Since  rf  coils  are  normally  the  complementary  part  of  a  tuned  circuit  whose 
function  in  the  electrical  network  may  be  either  frequency  controlled  and/or  matching,  it 
is  important  to  be  aware  of  how  the  coil  behaves  with  various  plasma  conditions.  This 
can  be  accomplished  by  examining  the  changes  in  coil  inductance  AL  and  coil  resistance 
AR  in  terms  of  x.  Thus,  if  Eq.  (30  is  substituted  into  Eqs.  (8)  and  (9),  then  the 
following  equations  for  AL  and  AR  as  a  function  of  x  are  obtained: 


and 


AL  _  L -L'  =  T _ _ 

L  L  x4  +  2x2  +  l\flx  +  1 


AR  _  R  R  =  xZ(l  +  y/Tx) 

R  R  x4  +  2x2  +  2\f2  x  +  1 


(31) 


(32) 


These  ratios  arc  plotted  in  Fig.  4.  From  this  figure,  it  can  be  seen  that  the  change  in  coil 
inductance  is  small  for  small  values  of  conductivity  and  increases  to  a  saturation  value 


Fig.  3  Experimental  and  Theoretical  Data  Fig.  4  Plot  of  AL/L  and  AR/QR  versus  x 
for  the  Figure  of  Merit  Qj  for  a  for  a  Solenoid  Immersed  in  a  Plasma 

Solenoid  Immersed  in  a  Plasma 
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k2L  as  the  conductivity  is  further  increased.  In  contrast,  the  change  in  coil  resistance 
rises  from  zero  for  zero  conductivity  to  a  maximum  value  {at  x  =  1.91)  and  then 
decreases  toward  zero  as  the  conductivity  is  further  increased.  From  this  analysis,  it  is 
observed  that  the  coil  is  most  sensitive  to  the  surrounding  plasma  when  x  is 
approximately  equal  to  two.  Thus,  for  a  given  range  of  plasma  conductivities,  an  rf  probe 
can  be  optimized  for  measurements  in  this  range  by  proper  selection  of  the  frequency  of 
excitation  and  the  probe  radius.  That  is,  x  is  a  function  of  w  and  a  (external  variables) 
and  o  (an  internal  variable). 

3.4  CALCULATION  OF  COIL  PARAMETERS  WHEN  THE  PLASMA  IS  IN 
THE  INTERIOR  OF  THE  SOLENOID 

The  coil-plasma  configuration  to  be  treated  in  this  section  is  illustrated  in  Fig.  5. 
In  this  case,  the  plasma  is  located  in  the  core  of  the  coil.2  Unlike  the  preceding  example, 

« -  l  - - 

0  O  O  O  O  O  0-0 - r 

VWVVYXX^AAAAAAjUAXlt  g  YJ  f 

b 


k  -  a/b 


o  o  o  o 


o  o  o 

Coil 


Fig.  6  Two-Dimensional  View  of  Solenoid  with  8  Plasma  Core 


the  electric  and  magnetic  fields  are  more  readily  obtained  from  an  exact  form  of 
Maxwell’s  equation  rather  than  by  introducing  an  approximate  expression  for  the  gradient 
operator.  Thus,  if  it  is  assumed  that  (1)  the  fields  vary  predominately  in  the  radial 
coordinate,  (2)  tiie  magnetic  and  electee  fields  are  of  the  approximate  forms, 
respectively,  Bc  =  i*Bs(p)expOc»;t]  and  fie  =  i^E^CpJexpLjoit],  and  (3)  Eg  is  zero  along  the 
z-axis,  then  the  magnetic  and  electric  fields  which  satisfy  Maxwell’s  equations  in  the 
region  0  <  p  <a  and  -£/2  <  z  <  £/2  are3 


Be  -  ij  B* 


— 

I  bar  ya  +  j  bei  ya 


(33) 


2  A3  though  Alt  form  Is  mod  sometime*  a  a  conductivity  trsuducer,  it  li  mo* t  often  found  in  conjunction 
with  if  plum  generation. 

3The  development  of  theta  solution*  from  Maxwell’*  equation*  it  prewited  in  Appendix  II. 
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and 


Er 


*  Y  Ba 
VSfV’ 


(berj  yp  +  bei  i  yp)  +  (be i ^  yp  -  beri  yp) 
ber  ya  +  j  bei  ye 


(34) 


where  Ba  is  the  value  of  the  magnetic  field  at  p  =  a,  7  =  s/upo  o,  and  her  and  bei  are 
Kelvin  functions  (Ref.  1 1 ). 

The  magnetic  field  in  Eq.  (33)  is  the  field  in  the  region  (0<p  <  a),  which  results 
from  both  the  eddy-currents  (is)  and  the  primary  current  (ip ).  That  is, 

-4  -♦  — ♦  . 

Bc  =  Bs  +  Bp  (35) 

where  5S  is  the  magnetic  field  due  only  to  the  eddy-currents  in  the  plasma  and  Bp  is  the 
magnetic  field  due  only  to  the  current  in  the  coil.4  It  is  important  to  be  aware  of  this 
when  calculating  the  figure  of  merit  Q8  since,  in  this  calculation,  only  the  fields 
associated  with  the  eddy-currents  can  be  used  in  calculating  energy  stored  and  energy 
losses.  The  magnetic  field  inside  a  long  solenoid  located  in  free  space  is  simply 

Bp  =  iz  Ba  exp  tjoit]  (36) 

Thus,  the  field  due  to  the  eddy-currents  alone  is  obtained  by  substituting  Eqs.  (33)  and 
(36)  into  Eq.  (35)  and  solving  for  Bs.  That  is, 


-*  A 

Bs  =  i* 


(bei  yp  —  ber  ya)  +  j  (bei  yp  -  bei  ya) 


ber  ya  +■  j  bei  ya 


exi 


ptjwtl 


(37) 


The  time-averaged  energy  stored  (Ws)  is 

In  2n 


Ws 


1  1/2  cH  •  /r  o  * 

-Tf  s  j  V- 

— 0/ 2  0  0  ' 


+  f 


Ec  •  Ec*1 


pdpd<f>dz 


(38) 


g  .  B  * 

However,  when  a  »  coe,  then  _ i-1 — —  »  <E-Ea*  ;  thus  Eq.  (38)  reduces  to 

Ho 


m 

nl  B  2  £  p  [(her  yp  a  ber  ya)2  +  (bei  yp  -  bei  ya)5]  dp 
2p0  [ber2  ya  +  bei2  ya] 


(39) 


or  after  integration 

co  nl  B2  a2 


Wa  = 


2^2  Ho 


(ber  x  +  bei  x)  berj  x  —  (ber  x  — bei  x)  beii  x  +  (ber2  x  +  bei2  x) 

V2 


x[ber2  x  +  bei2  x] 


(40) 


4lt  is  this  type  of  separation  of  fields  and  currents  in  which  the  development  of  a  transformer  model  for  a 
coil-plasma  system  becomes  apparent  from  a  field  viewpoint. 
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The  energy  loss  per  cycle  (Wd)  is 

ffd  “Q  T  f  J  J  oEo’E« “  ("S")  *«  Rb 


] 


(41) 


Substituting  for  Ec  from  Eq.  (34)  and  integrating  yield, 


/2»r\  /to  Jr  £  B  J  «2\  (ber  x  +  bei  x)  belj  x  +  !ber  x  —  be!  x)  bar^  x 
ffd  "  (S')  (77177  /  L  xtbar2  x  +  bei*  xl  “**. J 


Thug,  tiie  figure  of  merit  Q,  of  a  plasma  located  in  the  center  of  a  coil  ii 

(bar  x  +  bai  x)  berj  x  -  (bar  x  —  bai  x)  balx  *  +^5s;(ber2  x  +  bei2  x) 


Q. 


(bar  x  4-  bai  x)  bai  jt  +  (bar  x  -  bai  x)  berx  x 


(43) 


Equation  (43)  is  plotted  in 
Fig.  6  along  with  some 
experimental  data.  From 
these  data,  it  is  shown  that 
the  theory  agrees  quite 
well  with  experimental 
data,  which  is  indicative  of 
the  ease  and  reliability  of 
predicting  the  behavior  of 
solenoids  in  the  presence 
of  conductive  media.  Fur¬ 
ther,  this  dose  agreement 
between  the  experiments 
and  theory  (see  Figs.  3  and 
6)  indicates  that  the  as¬ 
sumptions  used  to  simplify 
and  obtain  closed-form  so¬ 
lutions  to  Maxwell’s  equa¬ 
tions  (e.g.,  infinitely  long 
coil)  can  be  used  with 
confidence  for  the  finite 
length  coil  Immersed  in  a 
conductive  medium. 


Fig.  6  Experimental  and  Theoretical  Data  for 
the  Figure  of  Merit  Q,  for  Solenoid 
with  a  Plasma  Core 
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SECTION  IV 

CONCLUDING  REMARKS 

Although  the  work  reported  herein  has  been  concerned  with  cylindricaliy 
symmetric  systems  and  with  conductors  with  uniform  conductivity,  it  should  not  be 
implied  that  this  approach  to  the  coupling  problem  is  limited  to  these  cases.  The 
transformer  model  is  equally  valid  for  uniform  and  nonuniform  conductors  and  also  for 
other  than  cylindricaliy  symmetric  systems.  If  the  system  is  more  complex  than  outlined 
here,  the  difficulty  will  be  in  obtaining  quantitative  information  since  Maxwell’s  equation 
in  this  case  will  be  much  more  difficult  to  solve. 

The  parameter  Qs  has  a  fundamental  definition  in  terms  of  energies  and  will  still 
be  a  useful  parameter  to  transform  from  circuit  to  field  phenomena,  even  in  the  more 
complex  systems.  For  example,  it  can  be  seen  from  Eq.  (14)  that  the  power  input  to  the 
plasma  is  maximum  when  Qg  =  1  regardless  of  the  state  of  the  plasma  and  the  coil.  Also 
from  Fig.  4,  it  can  be  observed  that  the  plasma  behaves  like  a  resistive  element  for  lower 
values  of  Q, ,  and  for  values  of  Qs  greater  than  or  equal  to  ten,  the  plasma  appears  to  be 
inductive.  Thus,  the  value  of  Qs  for  a  plasma  or  other  conductive  media  can  be  used  to 
describe  or  predict  the  behavior  of  the  plasma-coil  system,  and  since  Qs  is  a  function  of 
the  plasma  diameter,  frequency,  and  plasma  conductivity,  it  is  possible  to  adjust  these 
parameters  to  match  the  desired  test  conditions. 

It  is  recommended  that  Q,  be  adopted  as  a  correlating  parameter  in  all  rf 
experiments  (either  in  plasma  diagnostics  or  plasma  generation)  because  (1)  it  is  universal, 

(2)  it  encompasses  and  groups  together  the  essential  parameters  of  the  plasma  and  coil, 

(3)  it  is  nandimensional  and  thus  does  not  depend  on  magnetic  field  strengths,  currents, 
etc.,  and  (4)  it  is  readily  expressed  in  terms  of  either  circuit  ot  field  phenomena. 
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APPENDIX  I 

PARALLEL-TUNED  CIRCUIT  AND  MATCHING  PROBLEMS 


Generally,  the  rf-excited  coil  is  the  complementary  part  of  a  tuned  circuit.  This  is 
true  whether  the  coil  is  used  as  a  source  for  intense  electromagnetic  fields  required  for 
electrodeless  plasma  generation  or  as  a  plasma  diagnostic  tool  (e.g.,  conductivity 
transducer).  Since  the  coil  is  often  used  in  one  arm  of  a  resonant  circuit  and  because 
impedance  matching  between  the  tuned  circuit  and  the  electronics  is  often  a  critical  area 
on  which  proper  operation  of  the  system  depends,  then  it  becomes  especially  important 
to  study  this  particular  aspect  of  the  plasma-coil  relationship. 

This  appendix  is  concerned  with  the  resonant  system  in  which  the  coil  is  in 
parallel  with  a  capacitor  as  indicated  in  Fig.  1-1.  The  impedance  Z'  of  the 
capacitor-inductor  parallel  circuit  is 

Z'  _  +  j  jq'-  g)*l'c  (q'z  +  n/Qll 

(Q'2  +  1){<o2L'C)2  +  o'2  u  -2w2L'C) 


The  parallel  circuit  is  said  to  be  tuned  (i.e.,  imaginary  part  of  Z’  is  zero)  when 


l'2 


c  = 


<a2L'  \Q'2  +  1 


(1-2) 


The  real  or  resonant  value  of  Z'0  is  obtained  by  substituting  Eq.  (1-2)  into  Eq.  (1-1).  That 
is 


Z'  =  Q'<JL' 


Q'2  +  l 


.'2 


(1-3) 


If  Q'  and  L'  in  Eq.  (I-I)  are  expressed  in  terms  of  the  free-space  coil  parameters  L  and  Q 


L'  =  L 


R'  =  R 


1  - 


k2Qs 


1+Q 


s 


k2QQc 


1  + 


1+Q. 


o)L 

R' 


^k2Qs 


plasma 

Fig.  1-1  Parallel -Tuned  Circuit  with  the  Coil  in  the  Presence  of  a  Plasma 


coLI' 
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and  in  terms  of  the  parameters  k  and  Qf  by  using  the  transformer  model  (see  Fig.  1-1), 
then  Z'o  can  be  written  as 


where 


Z0'  =  Z0  f(k,  0B) 


Z0  =  QwL 


and 


f(k,  Q.) 


Q,2  (1  -  2k2  +  k*)  +  2k2Q,/Q  +  1 
Q.2 +kaQQ. +  1  _ 


a-4) 

(1-5) 

a-6) 


The  expression  for  Z'0  is  plotted  in  Fig.  1-2  for  various  values  of  the  coupling  coefficient 
k  and  for  a  free-space  coil  Q  of  100.  It  can  be  seen  from  examination  of  Eq.  (1-4)  and 
Fig.  1-2  that  the  impedance  of  the  parallel-tuned  circuit  at  resonance  contains  terms 
which  depend  on  the  physical  characteristics  of  the  coil  (e.g.,  free-space  inductance, 
free-space  Q,  and  the  "coupling  coefficient  k)  as  well  as  the  parameter  Q,  which  is 
functionally  related  to  the  physical  and  electrical  characteristics  of  the  plasma.  Since  the 
proper  operation  of  the  entire  system  (i.e.,  the  electronics-coil-plasma  complex)  depends 
quite  strongly  on  “matching  phenomena,1*  it  is  important  to  examine  this  aspect  in  some 
detail.  It  should  be  emphasized  too  that  this  electronics-coil-plasma  complex  is  a  coupled 
system  and  that  optimization  of  the  system  demands  that  all  aspects  of  matching  be 
considered  as  a  whole  rather  than  independently.  For  instance,  it  is  possible  that  an 
arrangement  which  might  constitute  an  optimum  match  between  the  coil  and  plasma 
would  create  an  unfavorable  matching  relationship  between  the  electronics  and  coil  which 
could  result  in  an  overall  decrease  in  performance  of  the  system. 


Qs - ► 

Fig.  1-2  Plot  of  Normalized  Impedance  Ratio  Z'q/Z„  versus  Q* 
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There  is  a  specific  area  in  which  this  aspect  of  the  problem  often  creates 
confusion  between  various  investigators  of  electrodeless  discharges.  This  can  be  illustrated 
by  considering  the  impedance  Z'0  and  the  power  input  to  the  plasma  from  a  system  like 
that  shown  in  Fig.  1-1.  The  input  to  the  plasma  is 


Pp! 


asma 


(1-7) 


From  this  expression,  it  is  clear  that  maximum  input  power  to  the  plasma  is  obtained 
when  k=  1  (actually  the  coupling  coefficient  must  always  be  less  than  unity)  and  when 
Qs  =  1 .  That  is, 


(P 


plasma 


d-8) 


Thus  it  appears  from  this  analysis  that  the  coil  should  be  closely  coupled  to  the  plasma 
(i.e.,  the  ratio  of  plasma  to  coil  diameter  should  approach  unity)  if  it  is  desirable  to 
maximize  the  power  input  to  the  plasma  for  a  given  current  Ip .  (In  practice,  a  mismatch 
generally  results  in  a  reduction  in  Ip.)  However,  when  k  and  Q,  approach  unity,  the 
impedance  Z'0  can  drop  from  its  free-space  value  Z0  by  as  much  as  two  orders  of 
magnitude  (see  Fig.  1-2).  This  can  have  very  serious  consequences  on  the  operation  of  the 
generator  which  supplies  the  rf  power  to  the  load,  and  it  may  happen  that  this  value  of 
impedance  is  much  too  low  for  coil-electronics  matching. 


Since  the  Class  C  oscillator  is  most  often  used  for  rf  power  generation,  it  is  not  a 
ample  matter  to  determine  the  most  favorable  value  of  the  load  Z’0  a  priori.  The 
efficiency  of  operation  is  closely  linked  with  the  value  of  the  load,  and  experimentally  it 
is  generally  observed  that  values  of  Z'D  ranging  from  1 000  to  1 0,000  ohms  are  required 
for  proper  matching. 

The  point  to  be  emphasized  here  is  that  the  overall  matching  problem  must  be 
viewed  in  order  to  arrive  at  the  proper  rf  system.  Sometimes  in  rf  plasma  discharges, 
system  performance  has  been  improved  by  decreasing  the  coupling  between  the  plasma 
and  the  coil.  In  cases  like  this,  the  explanation  is  that  the  impedance  match  between  the 
coil  and  the  generator  was  not  satisfactory,  and  by  decreasing  the  coupling  between  the 
coil  and  the  plasma,  the  impedance  of  the  tuned  circuit  was  increased,  which  created  a 
more  favorable  matching  condition  and  improved  system  performance. 
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APPENDIX  II 

SOLUTION  TO  MAXWELL'S  EQUATIONS  FOR 
A  CYLINDRICALLY  SYMMETRIC  SYSTEM 


The  coordinate  nomenclature  which  is  used  in  this  report  is  shown  in  Fig.  II-l. 


z 


Rg.  11-1  L^ia^ram  Illustrating  Cylindrical  Coordinate  (p,  <f>,  z)  and  Unit  Vecton 
(ip,  i(p,  iz)  with  Respect  to  CartBslan  Coordinates 


The  curl  of  a  vector  A  in  cylindrical  coordinates  is 


ft 


A  A  — 

pi<f>  iz 

d  d 

d<f>  dz 

Az  _ 


(IH) 
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Maxwell’s  equations  are 


VXE  =  - 

at 

(11-2) 

VXB  =  p07  +  Pap¬ 
al 

(H-3) 

O 

CO 

II 

o 

(H-4) 

V  ■  E  =  q/e 

(II-5) 

When,  (1)  the  electric  and  magnetic  fields  have  a  simple  sinusoidal  time  variation  (i.e,, 
ei“‘),  (2)  the  system  is  neutral  (i.e.,  q  =  0),  (3)  the  constitutive  equations  are  applicable 
(i.e.,  T=  all  and  D  =  elf),  and  (4)  the  conduction  current  is  much  greater  than  the 
displacement  current  (i.e.,  a»coe),  then  Maxwell’s  equations  reduce  to 


VXE  - 

(11-6) 

VXB  =  p0  a  E 

(H-7) 

<5 

DO 

II 

o 

(II-S) 

and 

V  ■  E  =  0 

(11*9) 

If  it  is  assumed  that  the  electric  and  magnetic  fields  are  of  the  form 

F  -  E0  (p)  e'wl 

(IMO) 

and 

B  =  tz  B,  (p)  ^ 

ai-11) 

then  the  partial  derivative  (3/3 p)  in  Eq.  (II-l)  can  be  replaced  by  the  total  derivative 
(d/d p).  Thus,  Eqs.  (II-6)  and  (II-7)  become,  respectively, 


- *■  A  TdE^ 

VXE  -  t-r 

+  j  =  *~iz  [jwBj 

(H-12) 

and 

1  »  a  dBx 

vxb  - d; 

A 

=  i<£  W  E0 

(H-13) 

or 

dB  z 

(11-14) 

dp 

Likewise 

(11—15) 

dE0  1 

d2B  , 

ip  dp2 
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An  equation  which  contains  only  the  magnetic  field  Bz  can  be  obtained  by  substituting 
Eqs.  (II-4)  and  (11-15)  into  Eq.  (11-12).  That  is. 


dSB, 


i  i  B«  .  _ 

+ - - - }^0o  Bz 

P  dp 


0 


01-16) 


This  equation  can  be  transformed  into  Bessel’s  equation  if  the  substitution 
z  =  pe-*i/4  y/  cc>ju0 o  is  made.  That  is,  Eq.  (11-36)  becomes 


d2B, 

iz1 


+  z 


dB, 

dz 


z2B, 


(11-17) 


The  general  solution  to  Eq.  (11-17)  it 


B2  =  AJ0(2)  +  C  Y 0(z) 


(11-18) 


where  J0(z)  and  Y0(z)  are  Bessel  functions  of  the  first  kind  and  second  kind, 
respectively,  both  of  zero  order.  For  a  system  in  which  Bz  must  be  finite  at  the  origin 
(i.e.,  p  =  0)  and  has  a  value  Bz  =  Ba  at  p  =  a,  then  Eq.  (11-18)  becomes 


Bz  -  Ba 


J a  y/coptBcT  p) 

Jo  VcuPo'ci  ») 


(n-19) 


The  solution  to  Maxwell’s  equation  for  the  electric  field  is  obtained  from  Eqs.  (11-14)  and 
(11-19).  That  is, 


—  — 


1 


dB, 

dp 


J  iCe“TT1/‘4r  P ) 

J0  a) 

-  _ 


01-20) 


Often  it  is  more  convenient  to  work  with  functions  with  real  arguments.  This  can  be 
done  in  the  case  of  Bz  and  E$  in  Eqs,  (11-19)  and  (11-20)  by  using  the  following  equation 
of  transformation: 

Jnle^74  v)  =  ^  [bwa  v  +  j  bei*  v]  (II“2I) 


where  n  is  the  order  and  v  is  the  argument.  Thus,  Eq.  (11-19)  becomes 

ber„(\/cj/I7^ p)  +  j  bei„  (\/ <up0cr  p) 


> _ _ 

6  I  ber0  {yfuip^a  a)  -I-  j  bei0  {\j <Ofl0a  a 


r 

r 


and  Eq.  (11-20)  becomes 


E^ 


CJ 


Bn 


bert  {\/ajp0(jp)  +  j  beij  {\/cop0C  p) 


ber0  {\Jap0o  a)  +  j  b«i0  (V&J/i0<7  a) 


(II-22) 


(11-23) 
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